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In mathematics, the classical Plateau problem consists of finding the surface of 
least area that spans a given rigid boundary curve. A physical realization of the 
problem is obtained by dipping a stiff wire frame of some given shape in soapy 
water and then removing it; the shape of the spanning soap film is a solution to 
the Plateau problem. But what happens if a soap film spans a loop of inextensible 
but flexible wire? We consider this simple query that couples Plateau's problem to 
Euler's Elastica: a special class of twist-free curves of given length that minimize 
their total squared curvature energy. The natural marriage of two of the oldest 
geometrical problems linking physics and mathematics leads to a quest for the 
shape of a minimal surface bounded by an elastic line: the Euler-Plateau problem. 
We use a combination of simple physical experiments with soap films that span soft 
filaments, scaling concepts, exact and asymptotic analysis combined with numerical 
simulations to explore some of the richness of the shapes that result. Our study 
raises questions of intrinsic interest in geometry and its natural links to a range of 
disciplines including materials science, polymer physics, architecture and even art. 
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1. Introduction 

Soap films are ethereal and ephemeral. Yet their beauty is not and has fascinated 
the scientist and the mathematician for more than two centuries, since Lagrange 
and Plateau first linked soap films and minimal surfaces: special surfaces of zero 
mean curvature that minimize the area enclosed by a given contour. The classical 
problem of finding the surface with least area that spans a given rigid boundary 
curve was first formulated by Lagrange (1760). However, the subject was lifted 
to a whole new level by Plateau in the mid nineteenth century using a series of 
beautiful experiments that showed that a physical realization of these objects nat- 
urally arises through a consideration of soap films (Plateau 1849). Since then, the 
subject has inspired mathematicians (Osserman 2002, Colding & Minicozzi 2007, 
Morgan 2008), scientists (Thomas et al. 1988, Kamien 2002), engineers and even 
artists (Ferguson) and spurred extraordinary developments in all these fields. Over 
the past two decades there has been a renewed appreciation for the subject as 
links to problems as disparate as materials science, string theory (Maldacena 1998) 
and other areas of high energy physics and cosmology (Penrose 1973) have become 
increasingly apparent. 

Almost without exception, the focus of problems over the last two centuries has 
always been on determining the shape of a soap film with a given rigid boundary. 
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Figure 1. A physical realization of the Euler-Plateau problem. (A)-(D) Examples of min- 
imal surfaces bounded by an elastic line obtained by dipping (and removing) a loop of 
fishing line (0.3 mm diameter) in soapy water. The interplay between the forces due to 
surface tension of the film and the elastic forces arising at the boundary can be quantified 
via the dimensionless number 7 = crL 3 /a with a the surface tension, L the length and a 
the bending rigidity of the boundary. For small values of 7, the soap film takes the form 
of (A) a planar disk, but upon increasing 7 it buckles in the form of (B) a saddle, (C) a 
twisted figure eight and (D) a two headed racket-like structure. 

However, since experiments with soap films are often carried out by dipping closed 
wire frames into soap and then pulling them out, there is a natural generalization of 
this boundary value problem that is suggested by the following question: what if the 
soap film is not bounded by a rigid contour, but has instead a soft boundary such 
as a flexible, inextensible wire? This simple question links at the hip the Plateau 
problem to another classical problem in geometry, that of the Elastica, a curve of 
given length that minimizes the total squared curvature, first formulated by Euler 
(1844) and since then object of continuous mathematical investigation (Bryant & 
Griffiths 1983, Langer & Singer 1984, Mumford 1993). The natural marriage of 
two of the oldest geometrical problems linking physics and mathematics leads to a 
quest for the shape of a minimal surface bounded by an elastic line. The competition 
between the surface tension of the film, that acts to minimize the area bounded by 
the loop, and the elasticity of the boundary, that tries to minimize its deformation, 
dictates the shape of this composite structure. Since small loops are stiff, they will 
remain flat and circular when spanned by a film. On increasing the size of the loop 
past some critical value, the forces due to surface tension can no longer be sustained 
by the boundary which buckles out of the plane even as the film takes on the shape 
of a saddle. Dimensional analysis suggests that the critical size of the loop when 
it just buckles scales as (a/a) 1 / 3 , where a = EI (E Young's modulus and / area 
moment of inertia of the elastic filament) is the bending stiffness of the elastic line 
and a is the surface tension of the soap film. 

In Fig. 1A-D, we show a gallery of shapes obtained via a simple experiment of 
dipping a loop of fishing-line into a solution of water and dish soap, and pulling it 
out. The resulting shapes are quite varied: short loops span a planar disk (Fig. 1A), 
intermediate loop lengths cause the film to spontaneously twist out of the plane 
(IB) and eventually form a planar figure-eight shape (Fig. 1C). At the center of 
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the eight (where the boundary crosses) the surface normal rotates by 180°, and the 
spanning minimal surface is helicoid-like. Increasing the length even further leads 
to a a two-headed racket-like structure (Fig. ID). 

To understand these shapes and the transitions between them, we start with 
a mathematical formulation of the problem using a variational principle, given in 
section 2. We then briefly describe some general geometrical and physical aspects 
of the resulting Euler-Lagrange equations before describing the first symmetry- 
breaking bifurcation from a planar circular disc to a planar elliptical disc in section 
3. In section 4 we describe numerical approximations to the secondary bifurcations 
from the planar disk to twisted strips that resemble our experiments. In section 
5, we characterize these using an asymptotic theory using a minimal trial shape 
that captures the symmetries of the solutions seen in our computations, and finally 
conclude in section 6 with a discussion of some open problems that our study points 
towards. 



2. Mathematical Formulation 

To quantify the equilibrium shapes of these softly restrained soap films, we use 
a minimal mathematical model that poses the resulting forms as a results of the 
minimization of the following functional representing the total energy of the system: 

F[M]=a [ dA + (f ds(aK 2 +/3). (2.1) 

J M JdM 

Here M is a mapping of a two-dimensional closed disk representing a film of surface 
tension a bounded by a flexible curve of circular cross-section with bending rigidity 
2a. The boundary dM is thus treated as an Elastica immersed in three-dimensional 
Euclidean space R 3 , and j3 is a Lagrange multiplier, analogous to line tension, that 
enforces inextensibility of the boundary: 

L = <£ ds (2.2) 

JdM 

The assumptions that underlie our minimal model are as follows: (i) the soap film 
itself has negligible thickness and is isotropic in its tangent plane, with an energy 
proportional to its area, (ii) the soap film meets the the bounding elastic filament 
at its edge along a line, so that it can apply normal forces to the filament, but 
no torques, (iii) the bounding elastic filament has a circular cross-section with a 
uniform bending and twisting stiffness along its entire length and its ends are glued 
together without twisting them relative to each other. 

Then, it follows that the filament must remain untwisted even when buckled 
and bent out of the plane, although it can and does have geometric torsion. This is 
because the translational symmetry of the circular cross-section along the filament 
implies that the twisting strain along it must be a constant (Love 1927); indeed the 
twist is the Noetherian invariant associated with this symmetry. Taken together 
with the fact that the closed filament has zero twist initially, this implies that 
the filament must always remain twist-free; therefore the energetic contributions 
associated with twisting deformations vanish identically. As we will show, our simple 
theory suffices to capture the qualitative nature of the shapes and the transitions 
between them. 
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To derive the conditions for the extremum of the functional (2.1), we will impose 
a small virtual deformation of the film and its boundary and use this to calculate 
the corresponding change in energy and the necessary conditions for the resulting 
shape to be a local minimum. For the boundary curve and changes therein, we use 
a parametrization in terms of the Frenet frame {t, n, 6} of the tangent, normal 
and binormal vector at each point of the curve (Kreyszig 1991). This orthogonal 
triad evolves in space according to the Frenet-Serret equations: 

i! = kti , n' = —Kt + rb . b' = rn . (2.3) 

where k and r are respectively the curvature and torsion of the boundary curve 
and (.)' = d(.)/ds in Eq. (2.3) where s is the arc-length. A linear variation of the 
boundary energy due to a small deformation can be easily calculated and leads to 
the following form: 

S f ds (aK 2 + (3) = f ds [ol(2k" + k 3 - 2t 2 k) - (3k}£ + a f ds (Ak't + 2kt')t] 

JdM JdM JdM 

(2.4) 

where £ and rj are respectively the displacement along the normal and binormal 
vectors (see Mumford 1993 for a pedagogical derivation). The linear variation of 
the area bounded by the curve can be expressed in a standard coordinate basis 
g i = diR (with R position vector of the film and d{ partial derivative with respect 
to the z-th coordinate) as (Lenz & Lipowsky 2000): 

s[ dA= I dA(V-u-2Hw) = (f dsv u-2[ dAHw. (2.5) 

J M J M JdM J M 

Here u and w are respectively the displacement along the tangent plane of the film 
(expressed in the basis and its normal direction A/", H is the mean curvature 
of the film and v = t x N is the outward directed tangent vector normal to the 
boundary curve. 

Compatibility demands that the bulk and boundary variations must be consis- 
tent with each other. To understand the conditions for this, we now need to express 
the first term on the right-hand side of Eq. (2.5) in the Frenet frame of the boundary 
curve. This can can be done by letting: 

N = costf n + sintf 6, (2.6) 

where $ is the angle between the surface normal to the film and the normal to the 
boundary curve and will be herein referred to as contact angle. Standard algebraic 
manipulations then yields: 

v • u = r] cos i!)b — £ sin $n . (2.7) 

Enforcing the condition that the energy variations (2.4) and (2.5) vanish while 
satisfying the compatibility condition (2.7) yields: 

k- sintf = 0, (2.8a) 
2a 

2k r + kt + — cos i? = , (2.8b) 
2a 



H = 0. 



(2.8c) 
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Eqs. (2.8) can be made dimensionless by rescaling the arc-lenght s, the curvature n 
and the torsion r by the length L of the boundary. This identifies unambiguously a 
single dimensionless parameter 7 = aL 3 /a representing the ratio between surface 
tension and bending rigidity. In sections 2 and 3 we will see how tuning 7 (for 
instance by changing the length of the boundary) triggers a sequence of buckling 
transitions, giving access to a variety of shapes. Determining the shape of the soap 
film then corresponds to solving (2.8) subject to periodicity of the boundary curve; 
here, we note that it is not sufficient to impose periodicity on the curvature and 
torsion, but we must impose it directly on the coordinate description of the curve. 

3. Some general considerations and bifurcation from the flat 

state 

Before we proceed to understand the nature of the solutions of the Euler-Lagrange 
equations derived in the previous section, we provide some general considerations of 
the nature of the equations. From a geometrical perspective, since the projection of 
the surface normal N on the boundary normal and binormal direction are associated 
with the normal and geodesic curvature of the boundary, namely: tvji — Kj COS d and 
Kg = ^sin^, an alternative way to write the equation satisfied by the boundary 
curve (2.8) in terms of the extrinsic and intrinsic curvature of the boundary leads 
to 

«" + ^« 3 - (r 2 ^f]n-^^={) 1 (3.1a) 

2k't + kt'+£-^=0. (3.1b) 
2a k 

From a physical perspective, is interesting to note that Eqs. (2.8) can be also 
obtained from the classical equilibrium equations of rods in presence of a body force 
directed along v. These require the total force and torque acting on a line element 
of the rod to vanish at equilibrium and yield (Landau & Lifshitz 1986): 

F' = -K , M f = Fxt. (3.2) 

The vectors F and M can be expressed in the Frenet frame of the bent rod as: 

F = (an 2 - ft) t + 2<W n + 2anr b , (3.3a) 

M = -2anb. (3.3b) 

We see that while torque balance is automatically satisfied, taking the external 
force K — and differentiating the expression for F readily yields Eqs. (2.8). 

We now turn to a number of properties of the system that can be inferred from 
Eqs. (2.8) and (3.1) without knowledge of the exact solution. Multiplying Eq. (3.1b) 
by k allows us to recast the resulting expression as (n 2 r)' = —(a/2o-)n ni which 
integrated along the length of the closed boundary yields the following integral 
formula for the normal curvature n n : 

(j)dsK n = 0, (3.4) 
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valid for all values of the physical parameters. 

Furthermore, it is possible to consider a special class of solutions of Eqs. (2.8) 
for which the contact angle $ is constant along the curve. Then Eq. (3.4) can be 
used to prove that $ must then be 90°. To understand this latter statement we first 
need to recall that the normal vector n is undefined at inflection points (i.e. where 
n = 0). From Eq. (2.7) it follows that, in order for $ to be defined everywhere 
along the boundary, the curvature of the boundary must be non-zero everywhere, 
i.e. ft > 0. If $ is constant, writing K n = n cos d in Eq. (3.4) yields: 



but since k is strictly positive, this last condition is possible if and only if cos$ = 0, 
thus $ — ±7r/2: i.e. n n = and the surface normal must be perpendicular to the 
curve normal everywhere if $ is to be constant. This imposes strong constraints on 
the nature of the solution in this case as we will now argue. 

Curves having zero normal curvature everywhere are called asymptotic curves 
and their tangent vectors asymptotic directions. Every curve that lies on the plane is 
clearly asymptotic. On a minimal surface, on the other hand, there are two asymp- 
totic curves passing through each point where the Gaussian curvature K is non-zero 
and these curves always intersect at 90°. If the minimal surface contains flat points 
(i.e. isolated points where K = 0), these are traversed by n > 2 asymptotic curves 
crossing at ir/n and forming, in their neighborhood, n valleys separated by ridges 
(Koch & Fisher 1990). These observations suggest that the only simply connected 
minimal disk bounded by a closed asymptotic curve lies in the plane; otherwise the 
asymptotic directions in the disk would form a vector field tangent to the boundary. 
In a simply connected disk this vector must then enclose a singularity or vortex, 
in the neighborhood of which the structure of the asymptotic directions would dis- 
agree with the previous classification. Hence, we conclude that the only solution 
having constant contact angle is given by the flat disk. 

However, even in this case, it is possible for the soap film to be non-circular 
when the parameter 7 is large enough. To see this, we note that a trivial solution 
of Eqs. (2.8) is a circular disk whose radius R satisfies the equation: 



for all parameter values, with f) enforcing the condition L = 2ttR with L the length 
of the boundary. For large enough surface tension, however, we might expect this 
configuration to become unstable and the disk to buckle into a more complex shape. 
To test this hypothesis we analyze the stability of the disk with respect to a small 
periodic displacement in the radial direction SR = pk sin k(j), where <\> is the polar 
angle and pk a small amplitude. Expanding the energy (2.1) to second order in pk 
yields: 




(3.5) 



crR 3 + /3R 2 -a = 0, 



(3.6) 




k 2 (3 a(2k 4 - 5fc 2 + 2) 
~R + R? 



(3.7) 



Here the Lagrange multiplier j3 is given by Eq. (3.6) with R = L/2tt, so that: 
P = [a — a (L/27r) 3 ]/(L/27r) 2 . Introducing the dimensionless parameter 7 = aL 3 /a 
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and replacing j3 in Eq. (3.7), it is easy to verify that the coefficient of the second 
order term becomes negative when: 

7 > 16tt 3 (/c 2 - 1). (3.8) 

Thus the first mode that goes unstable is that associated with an elliptical defor- 
mation with k = 2, corresponding to the critical value 7* = 487T 3 . This elastic in- 
stability was first discussed by Love (1927) in the case of a circular ring constrained 
to lie on the plane and subject to a uniform normal pressure, and has since been 
the subject of rigorous bifurcation analysis (Flaherty et al. 1972). As 7 is increased 
still further, the non-circular shape eventually twists out of the plane. To under- 
stand these transitions, we resort to a combination of numerical and asymptotic 
approximations. 

4. Numerical simulations of the soap film shapes 

To explore the variety of possible shapes resulting from the solution of Eqs. (2.8) 
and the transitions between them, we minimized a discrete analog of the total 
energy (2.1). The soap film is approximated as simplicial complex consisting of an 
unstructured triangular mesh. The internal edges of the triangles are treated as 
elastic springs of zero rest-length so that the total energy of the mesh is given by: 

F[M]=a (svWv + hY, M 2 - (4- 1 ) 

vedM eeM 

The first sum runs over the boundary vertices and (s v ) = (s v + s v -\)/2 is the 
average of the length of the two edges meeting at v. The curvature of the boundary 
is calculated as — I tfj ~^"v — 1 \/(s v ) with t v and t v -\ the tangent vectors at v. The 
second sum in Eq. (4.1) runs over all internal edges. If the triangles are equilateral, 
this yields a discrete approximation for the soap-film energy with the spring stiffness 
proportional to the surface tension, i.e. k = \/?>o- /A. The choice of minimizing the 
squared length of the edges, instead of the area of the triangles, is motived by 
numerical stability. Replacing the spring energy in (4.1) with the sum of the area 
of the triangles, has the effect of shortening the range of interaction between the 
boundary and the interior to the single strip of triangles at the boundary. Most 
standard local optimization algorithms would then attempt to reduce the area of 
these triangles to zero, thus suppressing the interaction between the boundary and 
the interior. In the continuum limit, the energy (4.1) approaches that given by (2.1) 
for original problem, and the sequence of shapes obtained with this method is in 
excellent agreement with our experimental observations shown in Fig. 1. 

Numerical minimization of Eq. (4.1) using a conjugate gradient method leads to 
a variety of shapes depending on the values of a, fc, the boundary length L as well 
as the initial shape of the domain. An interactive gallery of the shapes obtained 
based on our numerical simulations is available on-linef and we invite the reader 
to explore these shapes. A first set of simulations were run by starting with an 
unstructured triangular mesh of 474 vertices bounded by an elongated hexagonal 
perimeter; the mesh points are started out randomly displaced transverse to the 

f http : //www . seas . harvard . edu/sof tmat/Euler-Plateau- problem/ 
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plane. As the spring constant k was increased, corresponding to an increase in the 
surface tension, the system moved through a series of transitions shown in Fig. 2. 
Beyond a critical value of the surface tension (relative to the bending stiffness, as 
discussed earlier), the planar disk buckled into a two- fold mode consistent with our 
linear stability analysis. The geometry of this shape has been rigorously described 
(see for example Flaherty et al 1972, Arreaga et al. 2002, Vassilev et al. 2009) 
and is related with problem of the collapse of a infinitely long cylindrical pipe 
under a uniform pressure difference. For kL 3 /a « 643, the planar elliptical shape 
transitions to a twisted non-planar saddle-like shape, as experimentally observed. 
This is accompanied by an increase in the normal curvature of the boundary and 
there is a progressive twisting of the central "waist" of the surface. A further increase 
in k leads to an increase in the twist until the surface normal undergoes a full 180° 
rotation across the waist for kL 3 /a « 740. The surface then reverts to a planar 
configuration in the shape of the number eight (see Fig. 2E). 

Once the surface reaches the planar figure-eight conformation an increase in 
the surface tension does not produce further conformational changes until kL 3 /a ~ 
2000. Beyond this value, the surface again becomes non-planar, self-intersects at the 
pinch of the eight, and the two lateral lobes start bending toward each other. The 
curvature of the lobes increases with increasing k until they intersect to produce the 
complex shape shown Fig. 2H. The transition to this latter surface is analogous to 
the transition from the figure-eight to the two headed racket-like structure observed 
in our experimental realization of the problem. However, in the physical experiment, 
self-intersection cannot occur and adhesion favors the formation of a line of contact 
between the two lobes (see Fig. ID). 

Not surprisingly, given the nonlinear nature of the governing equations (2.8), 
the solutions are not expected to be unique, as our stability analysis of the flat 
disk has already indicated. In our numerical simulations, we also find that the final 
state is a strong function of the initial configuration of the system. If the initial 
configuration is not elongated, as in say a triangular mesh bounded by a hexagonal 
perimeter, which is closer to being circularly symmetric, we find other possible 
stable configurations. As in the previous case, for small values of the surface tension, 
the system rapidly relaxes into a flat circular disk. However, for kL 3 /a ~ 855 
the system transitions to a saddle-like configuration that is the classical Enneper 
minimal surface. Upon increasing the surface tension still further, the curvature of 
the boundary and film becomes larger still and the surfaces self-intersect leading to 
the beautiful strucure shown in Fig. 3D. 



5. Asymptotic analysis 

The numerical data show that the transition from the planar two-fold symmetric 
shape to the non-planar twisted eight is consistent with a supercritical pitchfork 
bifurcation, with (|ft n |) ~ (7 — 7*) 1 / 2 , where the angular brackets stand for an 
average along the boundary curve. Indeed Fig. 4 shows a best fit of the numerical 
data at the onset of the transition, the exponent obtained from the fit is 0.489 ± 
0.005. 

To understand the origin of the pitchfork bifurcation, it is helpful to introduce a 
specific approximate representation of the twisted soap film in terms of the following 
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Figure 2. Some solutions of the Euler-Plateau problem and a bifurcation diagram. Shapes 
obtained by minimizing the discrete approximation (4.1) of the continuum energy of a 
minimal surface bounded by an elastic line (2.1) as a function of the dimensionless number 
kL 3 /a. For small values of the parameter, the flat circular disk of radius L/2tt has the 
lowest energy. Upon increasing kL 3 /a the disk first buckles into flat elliptical shape (see 
text). For kL 3 /a > 643, this elliptical disk deforms into a three-dimensional saddle- like 
shape, becoming increasingly twisted with an increase in kL 3 /a > 643 until the surface 
eventually becomes flat again, but now forming a fully twisted eight-shaped conformation. 
This twisted flat shape is a local energetic minimum for 740 < kL 3 /a < 2000. For still 
larger kL 3 jot the twisted surface self-intersects at the waist of the figure-eight and the 
two lateral lobes bend toward each other, while for very large kL 3 / a the surface exhibits 
several self-intersections leading to the complex structure (H). Below the shapes we show 
the bifurcation diagram for the system, characterizing the shape using the absolute normal 
curvature of the boundary integrated along its length (a measure of the amplitude of the 
instability) as a function of the bifurcation parameter kL 3 /a, along with the location of 
the transitions described above. 
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(A) (B) (C) (D) 




Figure 3. Further equilibrium solutions of the Euler-Plateau problem. An initial configu- 
ration of a triangular mesh bounded by a hexagonal perimeter gives the classic Enneper 
minimal surface rather than the figure-eight shape shown in Fig. 1. 




kL 3 /a kL 3 /a 

Figure 4. Characterizing the onset of twist and non-planarity. The left panel shows the 
average absolute normal curvature of the boundary (| hv n | ) as a function of kL 3 /a in the 
neighborhood of the onset of the transition to non-planarity. The solid line on the left 
figure is obtained from a best fit of the data. The exponent resulting from the fit is 
0.489 ±0.005, and is consistent with the instability being of a supercritical pitchfork type. 
The right panel shows the average Gaussian curvature of the surface (K) as a function 
of kL 3 /a in the neighborhood of the onset of the transition to non-planarity. The linear 
behavior is also consistent with our approximate analysis which yields the result shown in 
the text. 



one-parameter family of surfaces: 

x = r (1 + t 2 ) cos (j) , y = r (1 — t 2 ) sin <j) , z = t (^J^J sm ^ > (5-1) 

where r G [0, R] and <p G [0, 2tt] are the usual plane-polar coordinates. Here t G 
[—1,1] is a parameter that characterizes the family of shapes. Eq. (5.1) is the 
simplest parameterization of a twisted saddle with disk topology and whose degree 
of non-planarity is controlled by the parameter t. For t = 0, the surface given by 
(5.1) is a disk of radius R; as t becomes positive (negative), the surface becomes a 
right-handed (left-handed) twisted figure-eight, and for t = ±1 the surface reduces 
to a flat twisted figure-eight bounded by a lemniscate of Gerono. In the following 
we will exploit the topological equivalence between this family of surfaces and the 
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actual twisted soap film to study the global properties of the surface at the onset 
of the transition, when \t\ <C 1. 

The first fundamental form for the surface (5.1) is given by: 



cjff — t I 2t 



lr 2 = t 4 + 2t 2 
while its Gaussian curvature is given by: 

K 



(^) 2 +cos20-(^) 2 cos# 

(s) s * nZ ^ — sin20 
(^) 2 -co S 20+(^) 2 cos# 



o(t 3 



(5.2a) 
(5.2b) 
(5.2c) 

(5.3) 



For the boundary curve, its arc- length and squared curvature are given respectively 
by: 

ds 2 = R 2 [(1 + t 2 ) 2 -2t 2 (cos20-cos40)] (5.4a) 
2 _ 1 + t 2 (10 - 6cos4(/>) -t 4 (2-6cos20-2cos6(/))+t 6 (lO-6cos4(/)) + t 8 



R 2 [(1 + t 2 ) 2 - 2t 2 (cos 2(f) - cos40)] 3 



(5.4b) 



Knowing the properties of the parametrized surface and its boundary in terms of 
its metric Eqs. (5.2) and curvature (5.4) we can approximate the total energy of 
the system given by (2.1) as: 



R 



aR 3 









-*H'?)] 







(5.5) 



Finally, the condition of inextensibility for the boundary curve is given by: 

L = j> ds = ttR (2 + 2t 2 - t A ) + o(t 6 ) . (5.6) 

Using Eq. (5.6) in Eq. (5.5) to eliminate R and then taking the derivative of the 
approximate energy Eq. (5.5) with respect to t yields the following equation of equi- 
librium for the amplitude parameter t that characterizes the shape of the boundary: 



t (96tt 3 - 7) + - jt 3 = , 
o 



(5.7) 



We recognize this as the normal form of a supercritical pitchfork bifurcation and 
see that it implies: 



f ~( 7 - 7 *) 



*\l/2 



7 



96tt 3 . 



(5.S 



We note that the critical value of the bifurcation parameter 7* is consistent with our 
exact linear stability analysis, but our approximate parametrization of the shape 
allows us to go beyond this via a weakly nonlinear analysis. 
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The normal and geodesic curvature of the boundary are given by: 
2t t 3 

K n = — — sin 20 + — (3 sin 20 — sin 40 + sin 60) + o(t 4 ) , 
1 t 2 (l + 3cos20 - 5 cos 40) , 4n 

«. = ^ + ^ + <**), 

which at leading order in £ leads to: 

<M>«8*~(7-7*) 1/2 , (5-10) 

and confirmed by our numerical simulations as shown in the left panel of Fig. 4. 

Having characterized the curvature of the elastic boundary curve, we now turn 
to the bounded minimal surface. The Gaussian curvature K of the surface in the 
neighborhood of the onset of the bifurcation to non-planarity is given by Eq. (5.3), 

jdsK g + J dA K = 2tt . (5.11) 

The integrated geodesic curvature can be easily calculated at leading order in t: 

j> dsK g ^ 2tt(1 + 2t 2 ) , 

from which, on using (5.11) and (5.10), we find that the integrated Gaussian cur- 
vature (K) = / dA K is : 

<K)~-( 7 -7*), (5.12) 

in excellent agreement with the numerical simulations shown in the bottom panel 
of Fie;. 4. 



6. Conclusions 

Through simple physical experiments using loops of flexible thread that are de- 
formed by a soap film that is bounded by the loops, we have demonstrated a 
natural link between two of the oldest problems in geometrical physics and the cal- 
culus of variations, bringing together Euler's elastica and Plateau's problem. The 
union of these leads to the question of the shape of minimal surfaces bounded by 
elastic lines, a new class of questions at the nexus of geometry and physics. We 
formulate a simple geometric variational principle for the shapes of these surfaces 
and their boundaries and a corresponding set of coupled system of partial differ- 
ential equations for the double free boundary problem. A combination of scaling 
concepts, asymptotic analysis and numerical simulation allows us to characterize 
the qualitative nature of the solutions that result. 

However, many questions remain. While a minimal physical realization of the 
Euler-Plateau problem is a kitchen-sink experiment, its ramifications are likely to 
go far beyond this specific manifestation, just as the humble soap film and the 
elastic filament have been relevant for the study of matter not just at the everyday 
scale, but also for systems that range from molecules to black holes (Thomas et 
al 1988, Kamien 2002, Penrose 1973), while inspiring art (Ferguson). In addition, 



Minimal surfaces bounded by elastic lines 



13 



both Euler's elastica and Plateau's problem have inspired a variety of mathematical 
developments in topology, geometry, analysis and beyond (Bryant & Griffiths 1983, 
Langer & Singer 1984, Osserman 2002, Colding & Minicozzi 2007, Morgan 2008). 
Perhaps their fruitful union might also do the same ? Some immediate mathematical 
questions include proofs of existence, regularity of solutions as well as measures of 
the non-uniqueness of solutions particularly for large values of the only parameter 
in the problem 7 or its discrete analog kL 3 /a. In terms of applications, the Euler- 
Plateau problem is likely lurking behind almost every physical manifestation of 
the Plateau problem, itself an idealization; time alone will uncover these. Further 
generalizations that introduce the effects of finite shearing and bending rigidity in 
the soap film, the torques from the soap film's Plateau borders on the filament, the 
effects of twist deformations of the filament due to anisotropy and inhomogeneity 
of its cross-section etc. are just some avenues that might be worth exploring; indeed 
some of these arise naturally in problems of vertebrate gut morphogenesis (Savin 
et al. 2011) and elsewhere. 
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